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Using  linear  theory,  expressions  for  lift,  drag  and  moments 
are  obtained  for  hydrofoil  operating  near  the  free  surface  and 
under  transverse  gravity.  Existing  results  are  easily  obtainable 
for  special  cases  when  gravity  is  neglected,  and/or  depth  below 
the  free  surface  approaches  infinity.  Numerical  calculations 
are  performed  and  results  are  shown  by  figures. 
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1.  Introduction 


Since  the  renewal  of  the  Interest  in  hydrodynamics  of  hydrofoils 
because  of  their  potential  advantages  when  used  on  small  craft  [i], 
variety  of  papers  have  appeared  dealing  with  different  aspects  of 
supercavitating  flow.  A complete  solution  to  the  supercavitating  flow 
is  the  analysis  of  the  flow  in  three  dimensions,  when  hydrofil,  is 
operating  near  the  surface  under  the  effect  of  transverse  gravity.  This 
problem  has  not  been  solved  In  its  generality  though  various  simplified 
versions  In  two  dimensions,  have  been  solved.  Most  of  the  work  has  been 
done  by  considering  two  dimensional  models  because  of  readily  availability 
of  technique  of  analytic  function  theory  for  such  problems . Tulin  W 
studied  the  hydrofoil  and  its  associated  cavity  as  a slender  body  and 
used  techniques  similar  to  Biot  [3),  Stewart  M and  others.  For  limited 
range  of  cavitation  numbers,  Tulin  used  this  theory  to  calculate  the  forces 
on  a flate  plate  hydrofoil  operating  at  small  angle  of  attack.  Wu  to  repre- 
sented an  extension  of  the  classical  free  stream  line  theory  of  Helmholtz 
and  Kirchoff  and  also[&  } extended  the  linearized  theory  to  include  the 
effects  of  camber  upon  the  forces  on  a profile  In  a cavity  flow.  These 
results  when  compared  with  experiments  £7  ] show  certain  small  but  symmetric 
discrepancies  for  the  circular  arc  and  the  flat  plate  profile  in  fully 
cavitating  flow.  It  may  be  noted  that  in  all  above  mentioned  studies, 
effect  of  gravity  and  free  surfaces  are  neglected.  Parkin  [ 8 ] studied 
the  cavity  flow  under  gravity,  possibility  to  account  for  the  discrepancies 
observed  in  Wu  [A  A)  theoretical  and  Parking J experimental  results. 

* Number  In  brackets  designate  references  which  are  listed  at  the  end 
of  this  work. 


Supercavitating  flows  under  gravity  has  also  been  studied  by  Acosta  j, 
Ldnau  and  Street  jjO  J , Street  M and  others.  Problem  involving  free  su- 
rface for  uniform  folw  is  treated  by  Auslaender  2J  and  for  non-uniform 
flow  by  Hsu  ["0  . Untill  now  problem  involving  gravity  and  free  sur- 
face has  not  been  discussed.  In  the  present  work  such  a problan  is 
undertaken. 


2.  Formulation  of  the  Problem 

Consider  an  inviscid  steady  two-dimensional  uniform  fluid  flow 
of  density  f . Suppose  the  flow  is  disturbed  by  introducing  a flat  plate 
of  unit  length,  inclined  at  an  angle  % to  the  direction  of  the  flow. 
(The  unit  length  body  is  used  without  loss  of  generality,  since  its 
use  is  equivalent  to  normalizing  the  problem  on  actual  length  of  the 
plate) . The  upstream  velocity  far  from  the  plate  is  U and  the  origin 
of  the  x,  y-coordinate  systan  is  taken  at  the  leading  edge  of  the  plate. 
The  supercavitating  hydrofoil  at  a finite  depth,  h , below  the  free 
surface  is  shown  in  figure  1.  It  is  assumed  that  : 

(a)  Hydrofoil  and  cavity  are  equivalent  to  a long  slender  body  which 
causes  only  small  pertubations  in  the  flow  that  is  uniform  far  from 
the  body. 

(b)  The  cavitation  number  is  zero 

Accordingly  the  velocity  components  are  assumed  to  be 

U (1+  u,v)  (1) 

where  the  components  of  the  perturbation  velocity  (u,v)  and  their 
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first  partial  derivatives  are  "small",  in  the  sense  that  their  square 
and  product  are  negligible.  These  components  may  be  combined  into  the 
total  complex  velocity  W(z)=  U(l+w(z)),  which  is  an  analytic  function 
of  z.  In  this  expression  the  complex  velocity  of  small  perturbation 
w(z)  is  given  by 

w(z)  = u - iv 

As  is  discussed  by  Auslaender  [lijand  3nith  the  simply  connected 

(Pig. 2)  polygon  A^  Bj*,  E F C**  can  be  mapped  into  the  lower  half 
of  the  plane  (Fig.  3)  by  the  conformal  transformation 

z = a A^-^-  - C 1 + "k  > J (2) 

where  - 71  < arg  z 0 

4 (^> 


and  aA7l  = h 


The  above  relations  have  been  obtained  by  using  the  constraints  that 
x = 0 and  x = 1 corresponds  to  0,  and  !j.  = 1 respectively. 

Also  it  has  been  shown  that 

when  h — » 0 , r a > ° 

1 A » I 


when  h 


ce> 


x - 10  — * 

t 

- (3) 

a » <so 

A > 2a. 

x + 10  y 

- (4) 

x-  10  » 

, *>° 

3.  Lift,  Drag  and  Moment' 
Coefficients 


By  normalizing  the  expression  for  lift  with  respect  to  1/2  f U 


the  lift  coefficient  is 

A?  ii 


cl“  j 


X fu‘ 


where  is  the  pressure  difference  between  the  face  and  back 

of  the  hydrofoil. 

By  using  equation  (2)  and  Bernoulli's  equation,  we  get 

Since  the  lower  surface  of  the  hydrofoil  is  a flat  place  (y0(x)=-  £ x) 
the  above  equation  reduces  to 


t «-i(  atv  i* 

y o J ^ 


+ 


(5) 


jl  $ + <*.  ' uz 

To  find  the  value  of  u(  \ ),  we  use  the  relation  between  v(  $ ) and 
u(^  ) as  discussed  by  Tulin  and  Burkart  [l^J,  we  have 

_>  1 f 1 u.  ( > , » , . . 

v(  T ) = Hr  1 1 — Z~, ® > where  0 ^ ^ ^ 1 

K 'o  V ' 


(b)  Angles  are  not  affected  by  the  transformation,  therefore 


ff-’W-T(t) 

(c)  For  the  flat  plate  = - £ 

Using  these  observations,  equation  (6)  become 


u (*) 


u ^ ,/x 


i*. 

o-V>*  <*-*') 


or 


u (5  ) = 


S ( 1-  2 


(7) 


substituting  value  of  u(^  ) as  given  by  the  equation  (7)  in  equation  (5) 
we  get 

CT 


■ + - o^]  + ^ 


By  using  equations  (4)  and  ( 8 ) give 


"L 


- Sn 

, 9s 

w - O 

< ^ 

f1]2.  ii 

] <K 

V. 

J u‘ 

TT  $ 

- 2 

1 

h=  <» 


(9) 


In  the  special  case,  when  the  effect  of  gravity  is  neglected, 
equation  (9)  gives  the  same  result  as  given  by  Auslaender  M- 
Also  for  infinity  depth  and  without  gravity  effect,  equation  (9) 
gives  the  same  result  as  given  by  Parkin J which  is  also  a linearized 
version  of  the  result  given  by  Lamb  M- 
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'fhe  drag  coefficient  can  be  obtained  from  the  consideration  that 
for  a fully  cavitating  flow  past  the  flat  plate,  the  local  pressure 
is  everywhere  normal  to  the  plate|s*J  • Therefore,  the  drag  coefficient 
( 0^  ) and  Cj  satisfy  the  relation  C / = tan  $ "=  % - — (10) 

using  equations  (9)  and  (10),  we  get 


K o + 


h = O 


U1 


CD  = 


l*  [p+  o'4-  £ ] + iiz  o -w  < oo 


-(11) 


Ti  S1  + 3 ^ 


ul 


W-  OO 


Moment  about  the  leading  edge  : 

By  dividing  the  usual  expression  for  moments  about  the  leading 
edge  by  1/2  PU  , we  get  the  coefficient  of  moment  . 


°M  = ~ 


pu1 


Ap.  X • J 


By  using  Bernoulli's  eaations  and  equations  (2)  and  (10) 


Si  3 ul  ~ 2 ^ ^ 


JL  J1 

a L 


V 

a + \ 
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By  carring  out  the  integration,  we  get 

A1  S 


CM 


a 


[-^-(2  0.-0  - a?  1 

l 2 Tt  J a +a  j J 

_2  U 

3 TJl 

- - (12) 


Equations  (8),  (11)  and  (12)  give  us  the  expressions  for  Cp,  CL  and 
. These  expressions  can  be  given  in  terms  of  the  Froude  number 
F = uAghpP  , based  on  the  hydrofoil  depth,  h,  as  follows  : 


CL  = 


sV. 


r,  a i4  i * 

|^C  »+»■>  - a-  J + 


Si 


_A J 

CL 


■ZAtl 


~ TT, 


!A“~  ) 1 - A 
A + A ' J 3 


4.  Numerical  Calculations  and  Results 

Before  making  any  numerical  calculations,  it  may  be  mentioned 

that  equations  (9)  and  (12)  conform  with  the  known  observation  that 

gravity  has  a significant  contribution,  when  velocity  is  small  and 

its  effect  may  be  neglected  for  large  velocity. 

Numerical  calculations  were  carried  out  on  IEM  370/125  and 

results  are  shown  in  figures  4 and  5«  The  effect  of  Froude  nurrber 

is  shown  in  figure  4 whore  the  variation  in  the  ratio  of  the  lift- 

coefficient  at  depth  h,  to  the  lift  coefficient  at  infinity  depth, 

i.e.,  CT/C  with  respect  to  F is  plotted  for  several  values  of 

h.  It  may  be  noted  that  for  hydrofoil  operating  at  a depth  greater 

than  1/2  of  the  chord  length,  the  value  of  CT/C,  remains  practically 

n l,  OO 

unchanged,  when  the  value  of  F decreases  from  large  values  to  F=  3. 
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This  result  is  similar  to  the  theoretical  results  of  Hough  and  Moran 
['*]  and  the  experimental  result  of  Strandhagen  and  Seikel  VS]  . 

The  main  difference  being  that  in  the  range  10 >F  >3,  which  corr- 
esponds to  the  most  foil  ooerating  conditions,  CT/Cr  does  not  decrease 

b 

as  rapidly  with  decreasing  F,  as  is  determined  by  Hough  and  Moran  [15). 
For  example,  when  h=  0.5  and  F decreases  from  10  to  3,  the  decrease  in 

value  of  C /C  at  F=3  is  only  2.9?  of  the  value  of  C /C  at  F=10. 

L b m,  L L oc 

Tt  may  also  be  seen  that  when  hydrofoil  is  operating  near  the  surface 
(h=0.1),  there  is  a sham  increase  of  34. 6?  in  the  value  of  C /CT 

L 1)  oo 

when  F increases  from  0.2  to  7 and  when  F increases  from  7 to  10, 
the  increase  in  value  of  C^/C,  is  less  than  1?.  Tt  can  also  be  safety 

OO 

concluded  that  for  all  practical  purposes,  the  value  of  CT/C  , for 

b L oe> 

large  value  of  F is  the  same  as  its  value  in  the  neighborhood  of  F=5, 
when  hydrofoil  is  operating  at  a depth  greater  than  1/2  of  the  chord 

length.  The  variation  of  C.-/C  with  h for  fixed  F can  be  easily  cal- 

b L eo 

culated.  It  is  clear  that  when  F increases  from  5 to  10  the  increase 
in  value  of  C /C  is  only  0.73? » 0.38?,  0.16?,  0.08?  respectivlty 

1 » J ' OO 

at  h=0.5,  h=l,  h=2  and  h=3  . This  again  confirms  that  for  all  practical 

purposes,  values  of  CT/CT  for  large  value  of  F is  the  same  as  at  F=5. 

Li  L'  0C> 

Tt  is  interesting  to  note  that  for  fixed  F (=5),  the  value  of  C /C 

b Lao 

decreases  steady  by  nearly  1.4?  for  each  Increment  of  0.1  in  value  of 
h,  as  h increases  steady  from  0.1  to  1.2  . There  is  a sharp  decline  in 
the  value  of  CT/C  of  5.5?  as  h increases  from  1.2  to  1.5,  when  h 

b La, 

increases  further,  value  of  CL/C  levels  off  with  its  value  close  to 

b b 00 

h=1.5  . Similar  pattern  follows  for  other  fixed  values  of  F . 
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The  results  for  the  variation  of  C / < with  respect  to  F,  for  fixed 

M * 

h are  shown  in  figure  5.  Figure  5 indicates  practically  no  change  in 

value  of  C^/5  as  F decreases  from  large  values  to  3 or  4,  a gradual 
decrease  in  value  of  C ^ & , as  F decreases  from  4 to  2,  and  a sharp 
decline  as  F decreases  beyond  2.  For  example,  when  h=0.97,  the  decrease 
in  value  of  Cy/g  is  2.8%  (of  its  value  at  F=  10)  when  F decreases  from 
10  to  4,  the  decrease  in  value  is  12.26%  (of  its  value  at  F=  4)  when  F 
decreases  from  4 to  2 and  the  decrease  in  value  is  42.9%  (of  the  value 
at  F=  2 ) when  F decreases  from  2 to  1. 
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Figure  5-  Variation  of  Cw / & with  respect  to  F,  for 

fixed  h . 


